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Abstract 

A quantity that promises to reveal important information on perturba- 
tive and non-perturbative QCD dynamics is the azimuthal decorrelation 
between jets in different hard processes. In order to access this information 
fixed-order NLO predictions need to be supplemented by resummation of 
logarithmic terms which are large in the region where the jets are nearly 
back-to-back in azimuth. In the present letter we carry out this resumma- 
tion to next-to-leading logarithmic accuracy explaining the important role 
played by the recombination scheme in general resummations for such jet 
observables. 



1 Introduction 



One of the most commonly measured jet observables in experimental QCD studies is the 
azimuthal decorrelation A0 between hard final-state jets. When compared to theoreti- 
cal estimates of the same, this quantity is expected to provide valuable information both 
on QCD parameters (strong coupling, pdfs) as well as dynamics in the near back-to- 
back region sensitive to multiple soft and/or collinear emissions and non-perturbative 
effects. To this end it has thus been often examined in experimental QCD studies at 
HERA and the Tevatron [1-4] , used for the tuning of parameters of Monte Carlo event 
generator models and to constrain unintegrated parton distribution functions (updfs) 
in conjunction with HERA data [5,6]. 

Various complementary approaches are possible to study the dijet A(f) including NLO 
calculations [7], resummation of logarithms arising from the back-to-back region A0 pa 
7r as well as non-perturbative effects important in the same region and in the small-rr 
regime the inclusion of BFKL effects or the use, as mentioned before, of unintegrated 
parton distribution functions [5,6]. However no attempt has been made at examining 
the use of a combination of these approaches in order to obtain general predictions valid 
over all of the dijet phase-space and widest possible range of A0. 

In ref. [3], for instance, a problem was noted in the comparison of NLO QCD cal- 
culations to HI data on the A0 distribution which was apparently resolved by the use 
of updfs [5] suggesting the importance of QCD dynamics and non-perturbative effects 
in the probed kinematical region. One is thus led to wonder whether other approaches 
based on say conventional resummation in the back-to-back region may also help to 
ameliorate the problems with pure NLO results. Such resummation would not include 
BFKL effects and it would be interesting to see at what x values genuine small- a; effects 
are actually needed by the data. 

However there has not been much progress in the case of final-state resummation for 
observables such as this which, in contrast to the much studied event-shape variables, 
are crucially dependent on the exact definition of final-state jets. These are typically 
quantities which are constructed from "aggregate" jet kinematic variables (momenta, 
azimuth, rapidities) of which other examples are dijet invariant masses and jet p t spec- 
tra. Here the jet momenta are obtained from the particle (hadron) momenta after 
running an algorithm and specifying a recombination scheme. The algorithm dictates 
which particles end up in the jet and the recombination scheme how the jet kinematic 
quantities are related to those of its constituent hadrons. Resummation in these cases 
is a far more delicate affair and there are only a few instances in the literature of 
resummed predictions for such jet observables [8,9]. 

One of the main complications that arises in such problems is that, as we shall il- 
lustrate, one is typically studying observables that are sensitive to energy flow outside 
well-defined jet regions which potentially means that many such observables fall into the 
category of non-global QCD observables [10,11]. Since it was shown that the resumma- 
tion of non-global observables is substantially more complicated than that for "global" 
quantities such as most event-shape variables and in any case restricted to the large 
N c approximation, the most accurate theoretical predictions can be obtained only for 
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global observables. This appears to rule out the possibility of complete next-to-leading 
logarithmic estimates for many interesting jet observables including potentially the az- 
imuthal decorrelation we study here. As far as existing predictions for jet observables 
are concerned, the issue of non-global logarithms was not dealt with in ref. [8] (pub- 
lished prior to the discovery of non-global effects) where they would arise in threshold 
resummation for one of the definitions (M 2 = (p 1 + P2) 2 ) of the dijet invariant mass 
studied there but would be absent for the definition M 2 = 2pi.p 2 . Further we should 
also mention here that the non-global component has been incorrectly treated in ref. [9] 
where it is mentioned that such effects will vanish with jet radius when in fact one 
obtains a saturation in the small R limit as was explicitly shown for the case of jets in 
ref. [12]. 

In the present letter we shall show an interplay between the potential non-global 
nature of the observable and the exact definition of the jet as provided by the choice of 
a recombination scheme. This may be taken as an example of how carefully selecting 
the definition of the observable and the jets one may be able to render an exact NLL 
resummation possible, avoiding altogether the non-global issue. To be precise, here we 
point out that in a certain experimentally popular recombination scheme (used to study 
dijet azimuthal decorrelations at HERA) the observable at hand is in fact global and 
hence one can resum up to next-to-leading logarithms exactly. In a different recombina- 
tion scheme (currently used at the Tevatron) the observable is non-global. However for 
the particular case of azimuthal decorrelations we point to recent developments which 
indicate that non-global logarithms while formally present in the latter scheme will be 
numerically insignificant here and should not substantially impede phenomenological 
investigations near the back-to-back region. We should mention explicitly that we do 
not advocate here the general use of one recombination scheme over another: a scheme 
that has good features theoretically for one observable may not be so good for another 
and hence ideally speaking an observable-by-observable choice is optimal. 

This letter is organised as follows. In the following section we derive the depen- 
dence of the quantity A0 on multiple soft emissions in two different recombination 
schemes and hence distinguish its global and non-global variants. In the subsequent 
section we provide resummed results in impact parameter space for the global variant 
for both DIS and hadron collisions while pointing out that we have also resummed 
the non-global variant to sufficient accuracy. We then present our numerical results 
for the inverse transform from impact parameter space and hence for the azimuthal 
decorrelation distribution. Lastly we briefly discuss our results mentioning the fur- 
ther developments needed in terms of matching to fixed-order calculations as well as 
including non-perturbative effects and point to work in progress in this regard. 

2 Recombination scheme, kinematics and global- 
ness 

We wish to study the impact of two recombination schemes used to construct the 
angle A0 between the final-state jets in dijet production. In the first scheme [13] 
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the jet azimuthal angle (pj is given by a pf-weighted sum over its hadronic constituents, 
4>j = Pt,i4>il Ylii^j Pt,ii while in the second scheme one constructs the jet four- vector 
Pj = YliejPii "with the sum running over hadrons in the jet, and then parameterises 
Pj = Pt,j {cosh rjj, cos 4>j, sin <f)j,smh.r]j) to obtain the jet azimuth (pj. The first scheme 
is employed for instance by the HI collaboration at HERA (see ref. [14]) while to our 
knowledge the latter (E'-scheme) is currently preferred by the Tevatron experiments. 

Having defined the relevant schemes let us consider the final-state kinematics. The 
final-state configuration that concerns us here is one where the hard jets are nearly 
back-to-back in azimuth and hence the system is close to the Born configuration for 
dijet production. In this limit other than the hard dijet system one has to consider 
the presence of any number of soft emitted quanta which cause a small deviation from 
— <pj2\ = A(p = tt. The transverse momenta of final-state particles can then be 
parameterised as belowS: 

Pt,i = Pt,i(l,0), 

Pt,2 = Pt, 2 (cos(vr - e), sin(vr - e)), 

= Pt,2(~ cose ' sine), 
h,i = k tji (cos (pi, sin 4>i), (1) 

where the hard final-state partons are labeled by 1 and 2 and the soft gluons by the 
label i. For only soft emissions the hard partons are nearly back-to-back and |e| <C 1. 

In the scheme involving the pt-weighted sum we write the azimuth of the leading 
jets as: 



Pt,i + Eieji Ki Pt 



<l>i* = ~ T V Z ^(vr-e) + J — , (2) 

Pt,2 + Z^j G j2 K t,i Pt 

where to obtain results correct to first order in the small quantities k t) i it suffices to set 
Pt,i = Pt,2 = Pt and by momentum conservation it follows that e = — k t ^ sin (pi/pt, 
discarding all correction terms quadratic in soft momenta, that do not affect our results. 

Note that the azimuth of the reconstructed jet 1 has a small deviation from = 0, 
whereas that for jet 2 has a small deviation from <p> = tt, due to the emission of soft 
gluons. Hence the effects of soft emission on the azimuthal angle (as measured by the 
deviation from A(p = tt) are given by: 



|tt - A0| 



k ■ 

} — (sin (pi - 6n(pi - 9 i2 (n - (pi)) 
^ Pt 



+ O (kt) , (3) 



where 9ij = 1 if particle % is clustered to jet j and is zero otherwise. The definition 
above implies that the observable in question is global since it is sensitive to soft emis- 
sions in the whole phase-space, both in and outside the jets, and the dependence on 



1 Here one is looking at the projections of particle momenta in the plane perpendicular to the beam 
direction in hadron collisions or that perpendicular to the 7* P axis in the DIS Breit or HCM frames. 
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soft emissions in either case is linear in k t . This property ensures that it is possible to 
resum the large- logarithms in the back-to-back region to next-to-leading (single) loga- 
rithmic accuracy without resorting to the large N c approximation needed for non-global 
observables [10,11]. 

Now turning to the i?-scheme, to obtain the corresponding dependence on soft 
emissions, we construct the four-momentum of a jet as p 1 ^ = J2 ie jPi, where the sum 
runs over all partons/hadrons in the jet. Thus one obtains, in particular, the transverse 
momentum vector of the jet and hence the angle (fij in the transverse plane as the 
inverse-tangent of the ratio of components, tan0j = pt, v /pt,x- Employing this procedure 
one obtains the following result for the azimuthal angle between jets in the four-vector 
recombination scheme: 

|vr - A(f)\ = 

where the sum extends over all soft particles not recombined with the hard jets. This 
result is natural since if all particles were combined into the hard jets then by momentum 
conservation the jets would be back-to-back in the plane transverse to the beam (no 
net transverse momentum). Hence with four- vector addition deviations from A0 = 7r 
are only caused by particle flow outside the jet regions. Observables sensitive to soft 
emissions in such delimited angular intervals are of the non-global variety [10,11], and 
hence in the i?-scheme definition of jets the azimuthal decorrelation is a non-global 
observable. 



— sin 0,- 

Pt 

j^jcts 



(4) 



3 Resummation 

Having established that the observable at hand is a global observable in the p t - weighted 
recombination scheme its resummation is now straightforward. It resembles closely 
resummation for the azimuthal decorrelations between hadrons studied in ref. [15] as 
well as the resummation of dijet rates in the region of symmetric E t cuts [12]. While 
refereing the reader to the above references for more detailed explanations of the steps 
involved, we sketch the main arguments briefly below. 

For the case of dijets produced in DIS one is examining soft radiation off a three- 
hard-parton antenna (taking account of the incoming parton in addition to the two 
hard partons that form final-state jets). One may thus describe the probability for n 
soft gluon emission by the essentially classical form [16]: 

\M n \ 2 = \M B \ 2 -f[W(k t ), (5) 

i=l 

with Mq the matrix element for the process with no emissions (Born order) and where 
one has the following gluon emission probability in the eikonal approximation: 

w (ki) = 9s Y {^i^ki) + w q2g (ki) - -^w qiq , 2 (h)*j , (6) 
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for emission of a soft gluon ki off a three-hard-parton system comprising quarks qi, 
qi and a gluon g. The Wy factors are just standard dipole antennae, with Wij(k) = 
(Pi-Pj) I {(jPi-k)(Pj-k)) representing the emission from the various dipoles formed by the 
three parton system, and g 2 = Ana s is the strong coupling^]. 

Given the factorised nature of the eikonal squared matrix element eq. (jSJ), it is 
only needed to factorise the phase-space in order to show exponentiation up to next- 
to-leading logarithmic accuracy. The phase-space constraint for computing the cross- 
section integrated up to some value of A0 can be expressed as the condition that one 
is considering events with | vr — A<p\ < A, where for small A one is in the back-to-back 
region. One can then obtain the A distribution or equivalently the distribution in A0 
by straightforward differentiation with respect to A. The integrated cross-section is 
(schematically) given by 



a(A)= J dB\M a B \ 2 /.(*, A*J) £ ^ / flWWaik 



a=q,g 



n ~ i=l 
n 



x A 



k 

) — (sin0i - Ba4>i ~ 0i2(it - (pi)) 



, (7) 



where we used a compact notation with the squared matrix element for lowest-order 
dijet production given by \M§\ 2 } and we integrate over the dijet configuration f dB in- 
cluding the experimental cuts, as well as denote by J[dki] the integration over soft gluon 
momentum components. Further the index a represents the type of incoming parton 
(quark or gluon) and f a (x,/j/i) the parton density with (if a factorisation sea 1<B In 
brackets we have the step function constraint that restricts real emission contributions 
while soft virtual emissions are unconstrained and will be included later by imposing 
unitarity. Since the emission probability for multiple soft gluon factorises as indicated 
above, to achieve a resummed result it only remains to factorise the phase-space con- 
dition by using a Fourier representation of the step function: 



A 



db 



sin 



(8) 



Using the factorised matrix element and phase space allows us to exponentiate the 
single gluon emission result in 6-space and one obtains: 



(t(A) = I dB\M£\ 2 Z a ({p},A) 



with 



Sa(M,A) 



a=q,g 

1 

7T 



db 



sm(6A) e -^/ a (2,/4 



(9) 



(10) 



2 The argument of the coupling is (to our next-to-leading logarithmic (NLL) accuracy) the transverse 
momentum of ki with respect to the dipole axis in the dipole rest frame k 2 ^- = 2/wij{k) [16]. 

3 x is not to be confused with Bjorken-x, it is the momentum fraction carried by the incoming rather 
than the struck parton and we have x = xb/^ where xb is Bjorken-a; and £ = Q 2 /(2p.q) with p the 
momentum of the incoming parton. 
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The function R a (b), known as the radiator, embodies the soft single-gluon result which 
exponentiates in 6-space. It contains a characteristic dependence on the hard parton 
configuration that we will presently explicate. This dependence is represented by the 
dependence on the set of Born momenta {p} of the function E above. 

Then we have for the radiator the standard-looking result: 

/d 3 k 
^-^-W a (k) (1 - eMMk)}) , (11) 

where we also introduced virtual corrections via the unity in parenthesis. Noting that 
to NLL accuracy one can replace (1 — exp[ibv(k))) — > 1 — cos(ibv(k)) — > 9 (y (k) — 1/6), 
where b = 6e 7B , we can write: 

£«(A,{p}) = -J - sin(6A)exp[- J R a (6)]/ a (x,fi}/P) , (12) 

and carry out the computation for i? a (6)0 

/d 3 k 
^y^-Wa(k) (v(k) - 6" 1 ) . (13) 

We have thus far accounted only for soft emissions. To extend the result to include 
hard collinear radiation we need to extend the computation of the radiator such that in 
the collinear limit we use the full QCD splitting functions instead of just the infrared 
pole pieces contained in the Wij antenna functions. Moreover a set of hard collinear 
emissions on the incoming leg are accommodated by a change of scale in the parton 
distributions f a (x, /^) — ► f a (x, l^ 2 f/b 2 ) via DGLAP evolution. Since this step is standard 
and common to 6-space resummations with incoming partons we do not display its 
derivation, but for a fuller treatment we point the reader to ref. [18]. 

The result for R a (b) including the extension for hard collinear radiation can be 
expressed in terms of three pieces each with a distinct physical origin: 

R a (b) = RUb) + R a out (b) ~ In S (5, M) , (14) 

with Rf a and R^ nt being the contributions generated by emissions collinear to the incom- 
ing (excluding the set of single-logarithms already resummed in the parton densities) 
and outgoing legs respectively. In addition to these jet functions we have a soft func- 
tion S(b, {p}) which resums soft emissions at large angles, and which depends on the 
geometry of the emitting hard ensemble expressed here as a dependence on the set of 
hard Born momenta {p}. 

While our results eventually include the two-loop running of the coupling which 
is necessary to obtain full NLL accuracy (compute the full functions gi and g-i), for 

4 i?(6) can be written in the well-known exponentiated form [17] R(b) = Lgi(a s L) + g2{<^ s L) + 
a s {g^a s L) + • • • , where L — In 6. NLL accuracy in our notation amounts to the complete computation 
of the g\ and <?2 functions. 
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brevity and to illustrate the main features we report our results here in a fixed coupling 
approximation. In this case we simply obtain: 



KutW = (C 1 a + C 2 a )^^ + -L^-ln3-41n2 + 31n^jj + 

4 n 

+ (%B«)L, (15) 
K(i>) = C?^(2Z, 2 + 4i(-ln2 + ln^) + 4C°^B»L, (16) 
ln 5 (Mrf) = -4L(2 C ^ln^ + C ^ln^), (17) 

with L = In 6. In the above Cf is the colour charge of the incoming parton in channel 
a, for instance Cf = Cp for a = q, the incoming quark channel. Likewise C" 2 are 
the colour charges of the partons initiating the outgoing jets 1 and 2 in channel a. 
The main aspect of the results for the collinear i?" ut in jet functions is a leading double 
logarithmic behaviour, where one notes the unfamiliar coefficient 2/3 (different from 
all commonly studied event-shape variables for instance) associated to the double logs 
on the outgoing legs, i.e. in the function -Rout- Additionally hard collinear radiation 
is described by single-logarithmic terms with the coefficients C^Bi for each leg, with 
the appropriate colour charge Cg (£ = i, 1, 2) and -Bj,i,2 depending on the identities 
(spins) of the incoming and outgoing partons such that Bi = —3/4 for fermions and 
B e = -{11C A - 4T R n f )/(12C A ) for a gluon. 

Finally we have the soft wide-angle single-logarithmic contribution In S, which de- 
pends on the geometry of the hard three-jet system via the dependence on dipole in- 
variant masses Qij = 2{pi.pj). This structure is characteristic of soft inter-jet radiation 
for three-jet systems (see e.g ref. [16] for a detailed discussion). The result can be easily 
extended to 2 — >• 2 hard processes as shown below. 



3.1 Radiator for hadron collisions 

One can easily generalise the results just presented to the case of azimuthal decorrela- 
tions in hadron collisions such as at the Tevatron and LHC Here one is dealing with 
the suppression of radiation from an ensemble of four hard partons since one has two 
incoming legs in addition to the final-state dijet system. 

One considers each partonic subprocess and obtains an equation similar to eq. (TT2"j) 
involving this time two pdfs for the incoming partons. The 6-space radiator can once 
again be computed dipole-by-dipole (i.e. for each pair of hard partons) and the results 
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combined after weighting by the colour factor for each dipole. It reads: 



Rout(b) 




Ql2 



Vt 



) 



) 



+ 



RUb) 



4 ry 

+-(C 1 B 1 + C 2 B 2 )^L, 
(C a + C a ) + AL (- In 2 + In ^) ) 



Vt J J 



+ 



(18) 



(19) 



In 5(6, M) 



+4 (Cji-Bjx + C i2 B i2 ) —-L, 
Tr (iJe-* rt / 2 Me-* r / 2 ) 



(20) 



Tr (HM) 



with £ = 2a s L/ir for a fixed coupling, and Ca and Cj2 being the colour charges for 
the incoming partons and C\ and C 2 those for the outgoing jets. The structure of the 
result is thus similar to that for the DIS (three hard parton) case except the function 
on the last line of the above equation, characteristic of soft wide-angle gluon radiation 
from an ensemble of four hard partons [19]. Here the quantity V is an anomalous 
dimension matrix while H consists of elements B^ representing the product of the 
Born amplitude in colour channel i and its complex conjugate in channel j. Lastly the 
matrix M represents a normalisation arising from the colour algebra. These matrices 
depend on the exact 2 — ^ 2 scattering channel considered as well as on the choice of 
colour basis (see e.g. ref. [19] for their explicit forms in particular bases). In the end one 
sums over all channels after folding the S(A) for each channel with the corresponding 
Born weights to obtain the final result. 

3.2 Non- global variants 

We have seen that in the .E-scheme, used at the Tevatron, the observable at hand is 
non-global and hence the resummation differs significantly from that detailed above. 
At the leading logarithmic level there are no double logarithms arising from the final- 
state jets, and so the collinear pieces proportional to the colour factors of outgoing jets 
in eq. (118]) would be absent. At the level of next-to-leading logarithmic terms, two 
additional pieces arise. One is the non-global piece computed for the two-jet case in 
for instance ref. [10], which is concerned with correlated multiple soft emission and can 
only be computed in the large N c approximation. The other piece is an "independent 
emission" contribution arising purely from the jet algorithm dependence which corrects 
eq. (ITS]) at the NLL level; in particular for the k t clustering algorithm this factor was 
seen to scale as -R 3 at leading order [20,21]. However, due to the fact that these effects 
arise first at O (a 2 s ), they become significant only in a region where the integrand in 
eq. (j!2p is numerically small [12]. For this reason they have a negligible impact on 
the resummation and can safely be ignored. Thus in practice we are able to provide a 
resummed result [22] also for the current experimental definition at the Tevatron. 
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4 Results and Discussion 



To provide a final resummed result for the A0 distribution one still needs to carry out 
the b integration in eq. ( fT2l) . In this section we describe how to produce numerically 
resummed differential distributions for A (and hence for A0) in DIS, starting from 
the resummed expression eq. (J9J), summed over the two incoming channels a = q,g (the 
generalisation to hadron collisions should then be obvious since the same considerations 
will apply there): 

*(A) =J2f dB \ M ®\ 2 s »( M> A ) • ( 21 ) 

«=<?,g 

We recall that the measure dB contains implicitly the acceptance cuts on the jet mo- 
menta, which in this case coincide with the Born momenta {p}. The function H a ({p}, A) 
is the NLL resummed distribution, which can be written as: 

S a (M, A) = - / x sin ( feA )/o e- R ^e- R ^S({p},b) . (22) 

7T J b 

At NLL level, the incoming and outgoing radiators i? in and R out and the soft function S 
are as given previously except that these are now re-computed with a two-loop running 
coupling, which is required to achieve complete NLL accuracy. 

However we now have to deal with an issue that has long plagued such resummations 
in 6-space in that some reasonable but perhaps somewhat ad-hoc prescriptions have to 
be adopted to practically evaluate the b integral in eq. (jUJ) (see for instance the energy- 
energy correlation in ref. [23] for related discussions): 

1. At large b the running coupling used to evaluate R(b) hits the Landau pole. 
We then decide to simply cut-off the b integral (i.e set R(b) = oo) at b = 
exp(l/(2a s (//R)/3o)) which corresponds to the Landau pole singularity. 

2. Additionally at large b, the factorisation scale of the parton density in eq. f[2"2"j) 
becomes small. Current parameterisations of parton densities usually fail for 
factorisation scales less than Qq = 1 GeV, corresponding to the fact that here it 
is not possible to neglect the intrinsic motion of partons inside the proton. This 
treatment is beyond the scope of the present letter, so we freeze the parton density 
at Q for b > fif/Q Q . 

3. At small b the radiator in eq. (iTTj) ought to vanish for 6 = 0, since this point 
corresponds to a complete cancelation between real and virtual correction terms. 
This cancelation is not present in the pure NLL approximation we obtain here. 
We have therefore set .Km/out (&) = S(b) = and freezed the parton density at /// 
for b < 1. Other modifications such as the replacement b — > Vl + b 2 that ensure 
a sensible 6 = behaviour can also be made but we have checked that this does 
not numerically alter our results presented here. 

An important thing to notice is that an NLL resummation is strictly valid for 1 < b < 
Hf/Qo. In particular the relative size of the large-6 part of the integral, corresponding to 
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120 130 140 150 160 170 180 

A <S? [deg] 

Figure 1: The resummed A(f> distribution for dijets in DIS. Also shown for comparison 
are the leading order (LO) and next-to-leading order (NLO) predictions from NLO- 
JET++ [7]. 



b > fif/Qo, can give us an idea of the impact of intrinsic parton transverse momentum, 
an area that we shall explore in more detail in forthcoming work [22] . 

We plot the resummed result for the A0 distribution in fig. [2] along with the fixed- 
order predictions for dijet production in DIS with Q 2 = 67 GeV 2 and xb = 2.86 x 10 -3 . 
These values and other cuts on the jets have been taken from the HI study [3] to 
which we would eventually compare our results. We fixed both renormalisation and 
factorisation scales to be the average transverse energy of the jets, and used CTEQ6M 
parton distributions [24], corresponding to a s {Mz) = 0.118. As we can see the fixed- 
order predictions diverge as expected near A0 = n. This divergence is cured by the 
resummation that goes to a fixed non-zero value at A0 = tt. Of note here is the absence 
of a Sudakov peak since the Sudakov mechanism does not dominate the b integral at 
very small A = \ir — A<f>\. The dominant mechanism to obtain back-to-back jets is thus 
a one-dimensional cancelation between emissions rather than a suppression of the k t 
of each individual emission, leading to a washout of the Sudakov peak as explained in 
detail in ref. [15]. 

In order to obtain complete predictions which can be compared to data two further 
developments need to be made. The first concerns matching to fixed order which is 
non-trivial since it requires information on the flavour of all the partons in the event 
which is not directly available in the fixed-order codes. Here we hope to exploit recent 
developments in this regard [25] which have addressed these issues in the context of 
hadron collider event shapes [26]. Secondly non-perturbative effects are expected to 
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play an important role in the region A0 w ir where they can be expected to significantly 
change the value of the distribution. Here one can apply a Gaussian smearing to our 
6-space results as is the practice for vector boson Q t spectra [27] as a model for non- 
perturbative effects whose parameters can be constrained phenomenologically. We leave 
both these developments for forthcoming work [22]. 
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